Astronomy & Astrophysics manuscript no. 8562 


February 2, 2008 


(DOI: will be inserted by hand later) 





o 
o 
o 

(N 

o 
O 



6 



SGR 1806-20 about two years after the giant flare: Suzaku, 
XMM-Newton and INTEGRAL observations 

p. s Mereghetti^, A. Tiengo^, S. Zane^ R. Turolla'^'^ D. Gotz^ N. Rea^-^, N. Kawai^ M. Ueno^ 

G.L. Israel^ L. Stella'', and M. Feroci^o 



Universita degli Studi di Pavia, Dipartimento di Fisica Nucleare e Teorica and INFN-Pavia, via Bassi 6, 1-27100 Pavia, Italy 

INAF - Istituto di Astrofisica Spaziale e Fisica Cosmica Milano, via Bassini 15, 1-20133 Milano, Italy 

Mullard Space Science Laboratory, University College London, Holmbury St. Mary, Dorking Surrey, RH5 6NT, UK 

Universita di Padova, Dipartimento di Fisica, via Marzolo 8, 1-35131 Padova, Italy 

CEA Saclay, DSM/DAPNIA/Service d'Astrophysique, F-91191 Gif-sur-Yvette, France 

SRON Netherlands Institute for Space Research, Sorbonnelaan 2, 3584 CA Utrecht, The Netherlands 

University of Amsterdam, Astronomical Institute "Anton Pannekoek", Kruislaan 403, 1098 SJ, Amsterdam, The Netherlands 

Department of Physics, Tokyo Institute of Technology, 2-12-1 Ookayama, Meguro-ku, Tokyo 152-8551, Japan 

INAF - Osservatorio Astronomico di Roma, via Frascati 33, 1-00040 Monteporzio Catone, Italy 

INAF - Istituto di Astrofisica Spaziale e Fisica Cosmica Roma, via del Fosso del Cavaliere 100, 1-00133 Rome, Italy 



Received / Accepted 



> 

00 

d 
o 



X 



Abstract In December 2004, the soft gamma-ray repeater SGR 1806-20 emitted the most powerful giant flare ever observed. 
This probably involved a large-scale rearrangement of the magnetosphere leading to observable variations in the properties of its 
X-ray emission. Here we present the results of the first Suzaku observation of SGR 1806-20, together with almost simultaneous 
observations with XMM-Newton and INTEGRAL. The source seems to have reached a state characterized by a flux close to the 
pre-flare level and by a relatively soft spectrum. Despite this, SGR 1806-20 remained quite active also after the giant flare, 
allowing us to study several short bursts observed by Suzaku in the 1-100 keV range. We discuss the broad-band spectral 
properties of SGR 1806-20, covering both persistent and bursting emission, in the context of the magnetar model, and consider 
its recent theoretical developments. 
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1. Introduction 

The four known Soft Gamma-ray Repeaters (SGRs) were dis- 
covered as transient sources of high-energy photons; they emit 
sporadic and short (~0.1 s) bursts of (relatively) soft gamma- 
rays with luminosity L ~ lO'^^'-lO"^^ erg s^'during periods 
of activity, that are often broken by long intervals of qui- 
escence. Three "giant" flares with luminosity >10"^^ergs ' 
have also been observed to date, each one from a different 
SGR: on March 5 1979 from SGR 0526-66 in the Large 
Magellanic Cloud jMazets et aljll979f) . on August 27, 1998 
from SGR 1900+14 ('Hurlev et al."l999f), and on December 27, 
2004 from SGR 1806-20 (Hurley et al. 2005). Persistent emis- 
sion with L ~ 10^^ erg s ' is also observed from SGRs in 
the soft X-ray range (<10keV) and, for S GR 1806-20 and 



SGR 1900-1-14, also in th e hard X-ray range dMereghetti et al 



2005al: iGotz et alj|2006h . In three cases, periodic pulsations 
at a few seconds have been detected. The bursts, the giant 
flares, the quiescent X-ray counterparts, and the pulsations 
have been interpreted in the framework of the magnetar model 



(see Thompson et alj|2002l and references therein). Magnetars 
are highly magnetized neutron stars with field strengths of 
lO'^^-lO'^ G, larger than those of the majority of radio pulsars. 
The ultimate source of energy for the bursts and the quiescent 
emission is believed to be the ultra-strong m agnetic field . 



SGR 1806-20 was discovered in 1979 dLaros etalJI 19861 



19871) and its persistent X-ray counterpart was observed for 



the fir st time with the ASCA satelHte in 1993 dMurakami et al 



19941) . A RossiXTE observation led to the discovery of pulsa- 
tions in the persistent emission with period P ^7.47 s and pe - 



2.6 X 10"^ s yr ' dKouveUotou et al.ll 19981) 



riod derivative P 
Under the assumption of pure magnetic dipole braking, these 
values imply a surface magnetic field strength of 8 x 10'^ G, 
strongly supporting the magnetar model. Both the burst rate 
and the X-ray persistent emission of SGR 18 06-20 started in- 
creasin g during 2003 and throughout 2004 dMereghetti et al 



creasin g durmg /UUj and throughout /UU4 (IMereghetti et al 
2005ci iTiengo et al.ll2005l IMereghetti et"al]|2007HWoods et al 



120071) . culminating with the giant flare of Dece mber 27, 2004 , 
during which -10"*^ erg were releasee^] CHurlev et al. 2005 



Send offprint requests to: Paolo Esposito, paoIoesp@iasf-miIano.inaf.it 



' Assumin g isotropi c luminosity and for a distance d = 15 kpc 
jCorbel et aU1997 : McCiure-Griffiths & Gaensleill2005h . 
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Mereghetti et al. I l2005bt iTerasawa et al.ll2005h . This giant flare 
was exceptionally intense, ~100 times more energetic than 
those from SGR 0526-66 and SGR 1900+14. Observations 
with RossiXTE unveiled, for the first time in an isolated neu- 
tron star, rapid quasi-periodic oscillations in the pulsating tail 
of the flare, likely relate d to global seismic oscillations on 
the neutron star surface ( Israel et al. 2005bf) . The flare pro- 
duced a hard X-ray (>80 keV) afterglow lasti ng a few hours 
jMereghetti et al.ll2005bl: iFrederiks ^ ^hOOi) and a radio af- 
terglow that faded in a few days (|Cameron et al. 2005). The 
small positional uncertainty of the radio observations per mitted 



to identify the likely IR counterpart of the SGR (iKosugi et al 
[2005; Israel et al. 2005a). The fluxes observed in the IR and 
gamma energy bands show a variabilit y correlated with that ob- 
served in the 2-10 keV energy range ( Mereghetti et alj|2007l) . 

After the giant flare, the persistent X-ray flux of 
SGR 1806-20 started to d ecrease from its outburst level, and its 
X-ray spectrum to s o ften jRea et al. 2005bla ;^ Mereghetti et al 

' J boovi 



2007t iTiengo et all l2005t IWoods et all f2007h . A Similar 



flux decrease have been observed from its radio afterglow 
( Gaensler et al. 2005; Taylor etal] 2005*) and its newly dis- 



covered IR counterpa rt (Ilsraeletalj r2005 a; R ea et alj|2005at 
Mereghetti eTal]|2007h . 



Here we present the results of the first Suzaku observa- 
tion of SGR 1806-20, covering both persistent and bursting 
emission in the 1-100 keV energy band. We also report on 
the analysis of a simultaneous observation performed with 
XMM-Newton and the latest outcomes of the monitoring of 
SGR 1806-20 with INTEGRAL, comparing them with what is 
seen in the same energy ranges with Suzaku. 

2. Suzaku observation and analysis 

The Suzaku observation of SGR 1806-20 started on September 
09, 2006 at 23: 13 UT and ended on September 1 1 , at 04:01 UT. 
The Suzaku X-ray observatory (Mitsuda et al ]l2007h carries on 
board the XIS spectrometers (Kovama et al .l2007h operatmg m 
the 0.2-12 keV energy ba nd, and the HXD collimated detec- 
tor (iTakahashi et al .l2007h . which covers the 10-70 keV energy 
range wifli PIN diodes and the 40-600 keV with GSO scin- 
tillators. Four X-ray telescop es with a spatial resolu tion (half- 
power diameter) of 2' (XRTs; ISerlemitsos et al ]E007 ) focus X- 
rays onto the four sensors (XIS 0, 1,2, and 3) that constitute the 
XIS instrument. Each XIS contains 1024 by 1024 pixel rows 
covering a 1 8' x 1 8' field of view, and features an energy resolu- 
tion of -140 eV at 6 keV. XIS 0, 2, and 3 are front-illuminated 
(FI) CCDs, while XIS 1 is a back-ifluminated (BI) CCD, that 
features an enhanced soft X-ray response. The XRT/XIS com- 
bination yields effective area per detector of ~330 cm^ (FI) or 
-370 cm^ (BI) at 1.5 keV, and of -160 cm^ (FI) or -1 10 cm^ 
(BI) at 8 keV. 

The 50 ks long observation was carried out with 
SGR 1806-20 at the "HXD nominal" pointing position. The 
XIS was operated in the normal mode with the 3x3 editing 
mode (time resolution of 8 s). The data sets were processed us- 
ing the version 6.1 of the FTOOLS package and the most recent 
available calibration files available at the time the reduction 
was performed (November 2006). The XIS pipeline products 



were affected by an imperfect charge transfer inefficiency (CTI) 
correction, resulting in a systematically lower energy scale. The 
error has been corrected by applying the CTI correction tool 
xispi again with correct CTI parameters. 

For the XISs, source spectra were extracted from cir- 
cular regions with radii of 3' centered at the position of 
SGR 1806-20, while the background spectra from composite 
regions (far enough from SGR 1806-20 to prevent contamina- 
tion by its photons). We screened the XIS events based on stan- 
dard criteri^ only events with GRADE 0, 2, 3, 4 and 6 were 
considered; the CLEANSIS script was used to remove hot or 
ffickering pixels; data collected within 256 s of passage through 
the South Atlantic Anomaly (SAA) were discarded; data were 
selected to be at more than 5° in elevation above the Earth rim 
(20° above the day Earth rim). This resulted in a net expo- 
sure time of 46.4 ks and about 57,000 net counts. The response 
matrices and effective area files were generated independently 
for each XIS with the tasks xisrmfgen and xissimarfgen 
(the ARE generator takes into account the level of hydrocarbon 
contamination on the optical blocking filter). The spectra were 
binned with GRPPHA following indications from the XIS Team; 
furthermore, the data were further rebinned to have at least 200 
source events per bin. 

The HXD data were selected according to the following 
standard criteria: at least 500 s after the SAA passages, day and 
night Earth elevation angles each >5°, and geomagnetic cut-off 
rigidity to be at least 8 GeV c"'. The exposure was corrected 
for the instrument dead time, for a net exposures of 48.4 ks 
in the PIN and 48.8 ks in the GSO. The HXD PIN and GSO 
instrumental background events were provided by the HXD 
Team (the instrumental background is due to events created 
by particles in the vicinity of the instrument). The source and 
background spectra (generated with the same good-time inter- 
vals) were both binned with GRPPHA following recommenda- 
tions from the HXD Team. 

To study the properties of the persistent emission of 
SGR 1806-20, we cleaned the event list from bursts by apply- 
ing intensity filters (with a negligible reduction of the net in- 
tegration time). Spectral fits were performed using the XSPEC 
version 12.3 software. The analysis of the bursts is presented in 
Sectionl23] 

2.1. Results in the 1-10 keV energy range 

Owing to the high interstellar absorption, very few counts 
were detected from SGR 1806-20 at low energies and thus we 
limited the spectral analysis to the 1.5-12 keV energy range. 
Suzaku is placed in a near-circular orbit around the Earth with 
an orbital period of about 96 minutes. Due to the source oc- 
cultation by the Earth in each orbit, the data-gathering required 
-1.2 days. Apart from the bursts, no variability in the XIS light 
curves of SGR 1806-20 was detected. The 8 s time resolution 
of the XIS data does not allow to detect the -7.6 s pulsations. 
We also investigated the possibility of spectral variability by 
splitting the observation in three intervals of equal duration, 



- See The Suzaku Data Reduction Guide, 
http : //suzaku . gsf c . nasa . gov/ docs/ suzaku/ analysis/ abc/ . 
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Table 1. Summary of the spectral results for the persistent emission in the 1.5-12 keV energy range with SuzakuPilS and 
XMM- Afewton/pn instruments. Errors are quoted at the 90% confidence level for a single parameter. 



ModeP 


Instruments 


Nh 
(10^- cm--) 


r 


(keV) 


(km) 


Flux"^ 
(10-" erg cm-- s-') 


(d.o.f.) 


PL 


XIS 


6.4 ± 0.2 


2.03 ± 0.04 






1.90 ±0.04 


1.16(283) 




pn 


6.5 ±0.2 


1.78 ±0.05 






1.74 ±0.04 


1.25 (71) 




XlS+pn 


6.5 ±0.1 


1.95 ±0.03 






1.73 ±0.03 


1.56 (356) 


PL + BB 


XIS 


7 1+0.6 
'■^-0.5 


1.8±0.1 


"•^^-0.07 


51 


2.1 ±0.1 


0.98 (281) 




pn 


"•'-0.7 




6+" - 

"•"-0.1 


2!] 


1 o+O.l 
^•^-0.2 


1.09 (69) 




XlS+pn 


6.9 ±0.4 


1.6 ±0.1 


0.55 ± 0.07 


T 7+1.6 
■-■^ -0.8 


1.8±0.1 


1.33 (354) 



^ Model applied in XSPEC notation: PL = phabs'-powerlaw and PL+BB = phabs*(powerlaw + bbodyrad). The abundances of 
[Anders & Grevesse (1989) are used throughout. 
Radius at infinity assuming a distance of 15 kpc. 
Flux in the 2-10 keV range, corrected for the absorption. 



with negative results. 

We fit simultaneously the XIS spectra (with relative nor- 
malization factors to account for the calibration uncertainties 
between the four cameras, see Section [3] for details) adopt- 
ing a power-law and a power-law plus blackbody model. The 
reduced of the former fit, x^r - 1-16 for 283 degrees of 
freedom (d.o.f.), corresponding to a null hypothesis probabil- 
ity of 0.03, is not completely satisfactory. The power-law plus 
blackbody model provided a better fit, with;^-^ = 0.98 for 281 
d.o.f. (null hypothesis probability - 0.6). The best fit param- 
eters are photon index F - 1.8 + 0.1, blackbody temperature 
keT = 0.49![J{J^ keV, and absorption A^h = 7.1+JJ-^ x lO^^ ^,^-2 
(Table [Til. The presence of the blackbody component is con- 
sistent with the fi ndings of deeper XMM-Newton observations 
dMereghetti et al.ll2007. see also Section|3]l. 



2.2. Results in the 10-100 keV energy range and 
broad-band spectral results 

The advantages of Suzaku/HXD over previous non imaging 
instruments are its small field of view (34'x34' FWHM be- 
low ~100 keV) and a low instrumental background. The im- 
ages obtained from INTEGRAL very deep exposures do not 
show contaminating point source s within the HXD field of 
view (see Mereghetti et al.l l2005ai Fig. 1) and no bright and 
hard X-ray sources below 10 k eV have been found either in 
the ASCA Galactic Plane survey ("Sug izaki et al. bOOl') or in the 
SIMBAD databas^E However, given that SGR 1806-20 lies at 
low Galactic latitude and longitude (b ^ 0° and I ^ 10°), the 
study of its emission in the hard X / soft gamma-ray band is 
complicated by the presence of the di ffuse emission from the 
Galactic Ridge (see lLebrun et al .|2004 and references therein). 

After standard data processing, a positive flux possibly as- 
sociated with SGR 1806-20 is detected in the HXD-PIN data 
up to ~40 keV (apart from the bursts, no significant emis- 
sion is detected in the GSO data). The instrumental back- 
ground counts obtained by simulations based on the present 
knowledge of HXD in-orbit performances, are about 70% of 
the ~26,400 total counts in the 12^0 keV band. To estimate 



the cosmic X-ray backgroun d level in the HXD-P IN band we 
took the spectrum reported in lGruber et al. (Il999h . of the form 

7.877 £-°-2'*e-^/41.13keV y.^y y.^y~l ^,^-2 ^-1 ^j.-!^ j^^^^jgi 

the Galactic Ridge emissio n we used the spectrum reported in 
VaUnia & Marshalll (Il998! ) for their Rl region (Central Ridge: 
-1.5° <b< 1.5° and -45° < / < 45°), where the SGR is lo- 
cated: a power-law with photon index F = 2.1 and surface 
brightness of 4.9 x lO"'^ erg cm"^ s"' 



sr-' in the 10-35 keV 



See http : // simbad . u- strasbg . f r/ simbad/ . 



band. 

However, the Rl region is wide and, since the Galactic 
Ridge emission str ongly varies with latitude and longitude 
(ILebrun et al.ll2004l) . this estimate of the Galactic Ridge con- 
tribution to the background could be severely inaccurate. 
Therefore, we analyzed a 43 ks long pointing carried out with 
Suzaku on April 07, 2007 at the coordinates b — Q° and 
Z - 8°, to provide a background field for t he observation of th e 
TeV source HESS J1804-216 (reported in Bamb a et al]|2007b . 
After a data processing and instrumental background subtrac- 
tion performed as described above, we fit the spectrum with 
a power-law with photon index fixed to 2.1. The ratio be- 
tween the measured power-law no rmalization and the one re- 
ported for the average spectrum in IValinia & Marshalll (Il998h 
is 0.9 + 0.3 (90% confidence level). So, these data indicate that 
the emission from the Galactic Ridge at a longitude closer to 
that of SGR 1806-20 is not significantly different from the av- 
erage value in the Rl region. We remark that this result relies 
on the accur acy of the instrument al background estimate, at 
present <5% ( Kokubun et al.l 2007 ). In fact, using a 5% higher 
background we obtain a nearly null Galactic Ridge emission, 
whereas a 5% smaller background yields a Galactic Ridge flu x 
2.0 + 0.4 times higher than that of IVaUnia & Marshalll (Il998l) . 

We note that to account for the whole signal detected in 
the HXD-PIN instrument, the Galactic Ridge emission in the 
HXD fiel d of view should be ~7 times higher than that re- 
ported bv' Valinia & Marshalll (Il998 '). This seems very unlikely 
to us and therefore we consider significant the detection of 
SGR 1806-20. However, both its spectral shape and flux are 
subject to the uncertainty reflecting the coarse knowledge of 
the Galactic Ridge contribution to the background. 

Including the cosmic diffuse and Galactic Ridge emis- 
sion as fixed components, we fitted the HXD-PIN spec- 



4 



P. Esposito et al.: X-ray observations of SGR 1806-20 



trum in the 12^0 keV band to a power-law model. 
The best-fit parameters are F = 2.0 + 0.2 and flux in 
flie 20-60 keV of (3.0 + 0.5) x lO"" erg cm^^ g-' (j^/^ = 
0.90 for 10 d.o.f.). The source, with a net count rate 
of ~0. 11 counts s"', accounts for ~70% of the counts re- 
mained after the subtraction of the instrumental background. 
Varying the assumed instrumental background by +5% we 
obtain best-fit fluxes of (3.7 + 0.5) x 10"" erg cm"^ s"' and 
(2.4 + 0.4) X 10"" erg cm^^ s"'. 

Fitting together the HXD-PIN and XI!^ spectra, we find 
that the HXD-PIN data must be scaled downward by a fac- 
tor of ~2 to be consistent with the parameters derived in the 
1.5-12 keV energy range. This scaling factor is unaccept- 
ably large, since the uncertainty in the relative calibration of 
the two instruments in the e nergy band considered here is of 
<20% dKokubun et al.ll2007h . To better reproduce the broad- 
band spectrum we tried a broken power-law plus blackbody 
model, with a normalization factor between the instruments 
kept at <1.2. We find an acceptable fit (x^ - 1.09 for 354 
d.o.f.; null hypothesis probability = 0.13) with the photon in- 
dex changing from 1.0 + 0.1 below the break at 16 + 2 keV 
to 2.2+°-4 above it, ksT = 0.8 + 0.1 keV, Rbb 
(at 15kpc), and A^h = 5.6+°;^ 
corresponding 2-10 keV and 20-60 keV unabsorbed fluxes are 
~2 X 10"" erg cm"- s"' and ~3 x 10"" erg cm"^ s"', respec- 
tively. The normalization factor assumed the value of 1.195, 
very close to the allowed maximum. 



2.5+°;4 km 



5.6+;!-j X 10^2 cm"2 (see Fig.[B. The 




E (keV) 



Figure 1. Broad-band Suzaku spectrum of SGR 1806-20 (see 
the online edition of the article for a color version of this fig- 
ure). The Suzaku'?, XIS 023 and HXD-PIN data (in black) are 
fit with the broken power-law (light blue dashed line) plus 
blackbody (red dot-dashed line) model. The XIS 1 data are not 
shown for clarity. We also plotted the INTEGRAL data (see 
Section|4]l using the blu circle marks. 



We did not find a significant pulsation in the HXD-PIN data 
(time resolution of 61 jUs). However, given the low signal-to- 
noise ratio, we do not expect to detect a clear signal if the 
pulsed fraction is ~10% (with a sinusoidal profile) as in the 
2-10 keV energy range, or smaller By folding the HXD-PIN 
light curve on the SGR 1806-20 pulsation period measured 
in simultaneous XMM-Newton data (Section |3]l and fitting it 
with a sinusoid, we determine a 3 cr confidence level upper 
limit of «20% on the amplitude of a sinusoidal modulation 
(Galactic Ridge emission subtracted). This upper limit is con- 
sistent with the preliminary results obtained with INTEGRAL 
in the 20-60 keV energy range (Gotz et al. 2007, in prepara- 
tion). 



2.3. Analysis of the bursts 

The Suzaku (XIS, PIN, and GSO) light curves show many short 
bursts (Fig.|2]l. To obtain significant constraints on spectral 
fit parameters, we considered only the events with more than 
50 counts in the HXD-PIN band. Thus, we selected only two 
bursts, indicated by the labels A and C in Fig. |2] and shown in 
Fig.S 



'j3 



C 
D 
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o 
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Time (s) 



0.8 
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Figure 3. HXD-PIN light curves of the bursts A and C in the 
energy range 12-50 keV. The origins of the time axis are arbi- 
trary and the time bin size is 10 ms. 



For the broad-band analysis we added the spectra of the three FI 
CCDs (XIS 0, 2, and 3) using the FTOOL matpha. We generated the 
instrumental responses by summing the redistribution matrices and the 
etfective area files using the FTOOLS addrmf and addarf . 



Since the frame-time of the XIS instruments (8 s) is much 
larger than the burst duration, we measured the duration of 
the bursts in the high time resolution HXD-PIN light curve 
(~0.5 s for burst A and ~0.4 s for burst C) and we set the 
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Figure 2. Temporal coverage of the Suzaku and XMM-Newton data (source plus background counts). Several bursts have been 
detected in the various instruments, but the brigh burst C is not visible in XMM-Newton/pn because it is too strong and saturated 
the detector. Despite this, the burst rate in the XMM-Newton/pn is higher than in in the Suzaku/XlS, with 0.5 + 0.2 burst ks"' 
versus 0.15 + 0.06 burst ks The difference is likely due to faint bursts undetected in the XIS owing to the long frame-time (8 
s). The Suzaku arrival times of the photons were barycentered to the Solar System using the task aebarycen. The comparison 
between the Suzaku/HXD and XMM-NewtonfEPlC times of bursts B and C proves a relative timing accuracy better than 0.2 s. 



Table 2. Burst spectral parameters for the two-blackbody model (assuming for the absorption the value of 6.5 x 10 
cm"^). Errors are quoted at the 90% confidence level for a single parameter. 



,22 



Burst 


ModeP 




^BB l'' 








r 


Flux"^ 


Xi (d.o.f.) 






(keV) 


(km) 


(keV) 


(km) 


(MeV) 




(erg cm"- s"') 




A 


BB +BB 


9 O+0.6 




10 ±2 


-'•-'-0.6 






1.3 X 10"' 


1.15 (47) 




BB +BR 


4 6+1.2 


5±2 








1.2 X 10-' 


1.08 (47) 




BB +PL 




-"-1 






1 5+"- 


1.1 X 10^' 


1.19(47) 


C 


BB +BB 


9 9+0.5 


20!5 


7 ± 1 


4± 1 




2.1 X 10^' 


0.95 (63) 




BB +BR 


A 7+0.5 


9±2 










2.0 X 10-' 


1.09 (63) 




BB +PL 


4.2 ± 0.4 


9±2 








1.4 ±0.2 


2.0 X 10"' 


1.12(63) 



" Model applied in XSPEC notation: BB + BB = phabs* (bbodyrad 
and BB + PL = phabs*(bbodyrad + powerlaw). 
Radius at infinity assuming a distance of 15 kpc. 
Flux in the 2-100 keV range, corrected for the absorption. 



bbodyrad), BB + BR = phabs* (bbodyrad + bremss). 



integration time of the XIS spectra to these valuefl The threshold above which the photon pileup becomes significant 

in the XIS cameras is ~100 counts frame"' for a point source, 

5 During the two 8 s frames containing bursts A and C, less than and both bursts had a count rate above this value. After a 

1.5% of the total XIS counts can be attributed to the persistent source ^.^gf^j analysis of the distribution of the counts in the pixels, 

plus background emission. Thus, the efi-ect of this emission on the XIS excluded a circular region with radius of 15" at the image 
spectra of the bursts is negligible. 
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center from the event extraction region to minimize the pileup 
effects. To increase the statistics in the spectra we added the 
four XIS spectra using the FTOOL matpha. We calculated 
the corresponding instrumental responses by summing the 
redistribution matrices and the new effective area files for the 
annular regions using the FTOOLS addrmf and addarf. XIS 
and HXD spectra were rebinned to have a minimum number of 
20 counts in each bin. For the background subtraction, under 
the usual working hypothesis that the bursting emission is 
present on top of the quiescent one, we used the spectra of the 
whole observation, cleaned from all the bursts. 

The column density of neutral absorbing gas along the line 
of sight is a critical parameter for the spectral fitting. The most 
precise measurement of this parameter comes from modeling 
the persistent X-ray emission: for SGR 1806-20 the column 
density is consis tently measured at A^h - 6.5 x 10^^ cm"^ 
(see Table[T] and iMereghetti et all 120071) . To test different 
spectral models we therefore decided to fix the equivalent 
hydrogen column at the value measured for the persistent 
source before and after the bursts. Spectral fits with power-law 
or thermal bremsstrahlung models, with the A^h value fixed at 
6.5 X 10^^ cm"^, yield unacceptable values. Good fits with 
either the power-law or bremsstrahlung, if the A^h is left free 
to vary, require large absorption: A^e values of 1.7 x 10^^ cm"^ 
and 1.4 X 10^^ cm~^, respective ly. Simil ar results were ob - 



tained by Fenimore et al. ( 1994t) and by iFeroci et aL (2004) 
for the soft gamma-ray repeater SGR1900H-14. A single 
blackbody spectrum gives a formally acceptable fit in the 
soft range, but it severely underestimates the observed flux 
at higher photon energies. Another possi ble thermal model 
is t he sum of two black bodies, as used by lOlive et al.l (120041) 



and Feroci et al 



(120041) to fit bursts from SGR 1900-H14. We 
found that this two components model provided good fits to 
the bursts, with the parameters summarized in Table |2] We 
tried other types of spectral models, obtaining almost equally 
good fits with either a blackbody plus power-law model or a 
blackbody plus thermal bremsstrahlung model (see Table |2] 
and Fig.lU. 



3. XMM-Newton observation and comparison with 
Suzaku 

The XMM-Newton observation of SGR 1806-20 started on 
September 10, 2006 at 10:11 UT and ended at 19:04 UT. It 
is therefore simultaneous to part of the Suzaku observation, as 
shown in Fig.|2] Here we present the analysis of the data col- 
lected with the EPI C instrument, which consists of t wo MOS 
(ITurner et al ] l200lh and one pn (iStriider et al.ll2001h cameras 
sensitive to photons with energies between 0.1 and 15 keV. 
The pn was operated in Small Window mode (time resolution 
6 ms), while the MOS 1 unit was in Timing mode (time resolu- 
tion 1 .5 ms) and the MOS2 in Full Frame mode (time resolu- 
tion 2.6 s); all the detectors mounted the medium thickness fil- 
ter. The data reduction was perfo rmed following the procedure 
described in lTiengo et alJ (2005) but using the XMM-Newton 
Science Analysis Software (SAS) version 7.0. 

By an inspection of the XMM-Newton lightcurves we found 



several bursts (see Fig.|2]i, but the XMM-Newton data do not 
provide significant improvement on the results obtained from 
the Suzaku data (Section [2.3b . In fact burst A occurred before 
the start of the XMM-Newton exposure, while burst C was too 
bright for the pn (the two MOSs registered few counts only), 
saturating the instrument telemetrjO. Only burst B was ob- 
served by both satellites (see Fig.|2]), but it was too faint for 
a spectral study. To obtain the results presented in this section, 
we excluded the bursts from the analysis by applying intensity 
filters. A further cleaning was necessary because of the pres- 
ence of soft proton flares during the observation. On the whole, 
the net exposure time was reduced from ~22.3 ks to ~21.6 ks 
for the pn detector. 

We fit in the 1.5-12 keV range the spectrum obtained 
with the pn camera adopting the same models used for the 
Suzaku spectral analysis. We again found that the power- 
law plus blackbody model provides a slightly better fit (see 
Table [U; the best-fit parameters are F = 1.6, kgT = 0.6 keV, 
and A^H = 6.7 x 10^^ cm"^ with a of 1.09 for 69 d.o.f.. 
The corresponding luminosity in the 2-10 keV band is 
5 X lO^^t/jg erg s"', where we indicate with d^^ the distance in 
units of N kpc. This luminosity is slightly higher than that of 
the pr evious XMM-Newton ob servation performed on April 4, 
2006 (IMereghetti et al.ll2007h . as also supported by a simple 
comparison of the pn net count rates (0.993 + 0.007 counts s ' 
with respect to 0.946 + 0.007 counts s"'), while the other spec- 
tral parameters are consistent. We also analyzed in a similar 
way the spectra obtained with the MOS cameras, finding re- 
sults consistent with the pn ones. 

Given the simultaneity of the XMM-Newton and Suzaku ob- 
servations, we have tried to fit together the XMM-Newton/pn 
and Suzaku/XlS spectra, with a normalization factor to account 
for the uncertainty in the absolute flux estimate of the differ- 
ent instruments. In Table [T] we report the results of such anal- 
ysis using either a power-law or a power-law plus blackbody 
model. We note that the values are unacceptable. Since no 
intrinsic spectral variability during the non-coincident expo- 
sure windows of the two satellites is expected (the bursts have 
been removed in both datasets), this simultaneous fit can be 
used to evaluate the cross-calibration discrepancies between 
the XMM-Newton/pn and the four Suzaku/XlS detectors. Fixing 
the XMM-Newton/pn normalization factor to 1 and linking all 
the other parameters, we derive the normalization factors re- 
ported in Table |3] either for the whole energy range or restrict- 
ing the fit to three energy bands. These values show that the 
four XIS detectors measure a systematically higher flux with 
respect to the pn, especially in the soft energy range. However 
we note that the cross-calibration accuracy between the pn and 
the XIS is, especially above 4 keV, of the same order of the dis- 
crepancies between the different XIS units. 

To derive the period of SGR 1806-20 we used the pn 
data. Photon arrival times were converted to the Solar System 
barycenter using the SAS task barycen. With a standard fold- 
ing analysis of the light curves, we measured a spin pe- 



* If the count rate in the pn is higher than the telemetry limit 
(~600 counts s"' for the imaging modes), then the so-called "count- 
ing mode" is triggered and for some time the science data are lost. 
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Figure 4. Suzaku broad-band spectra and residuals of the bursts A and C (XIS 1 .5-12 keV, HXD-PIN 12-50 keV and HXD-GSO 
50-100 keV). The model adopted is indicated on each panel (see also Table|2]l 



riod of 7.5891 + 0.0002 s. The resulting peak-to-trough pulsed Fmin are the observed background-subtracted count rates at the 
fractiorQ f Fpt = (Fmax - ^mm)/(^^max + ^^min), whcrc F^ax and peak and at the minimum, is (11 + 2)% in the 2-lOkeV band. 

^ No te that here we use a dif ferent definition of the pulsed fraction 
than in iMereghetti et all ( |2007|) . Using the old definition based on a 
sinusoidal fit to the profile, the estimated pulsed fraction in the same 
energy range is (8 ± 1)%. 
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Table 3. Normalization factors for the Suzaku/XlS cameras 
with respect to the XMM-Newton/pn (see Section[3]for details). 



are consistent with an extrapolation of the low-energy model 
(Fig.S. 



Detector 






Energy range 








1.5- 


12keV 


1.5-4 keV 


4-8 


keV 


8-12 keV 


XISO 


1.11 


±0.02 


1.23 ±0.04 


1.05 


±0.03 


1.1 ±0.1 


XISl 


1.08 


±0.02 


1.19 ±0.04 


1.00 


±0.03 


1.1 ±0.2 


XIS2 


1.17 


±0.02 


1.26 ±0.04 


1.13 


±0.03 


1.1 ±0.1 


XIS3 


1.19 


±0.02 


1.33 ±0.04 


1.12 


±0.03 


1.1 ±0.1 



4. INTEGRAL observations and broad-band 
spectral analysis 

We analyzed the INTEGRAL'S AO-4 Key Project (KP) ob- 
servation of the Galactic Centre, and report here results ob- 
tained with ISGR I, the low energy detector plane (15 keV- 

imager IBldl 



lMeV:LebrunetalJ 12003) of INTEGRALS 



dUbertini et al.ll2003h . The IBIS/ISGRI data were reduced us- 



ing the Off-line Scientific Analysis package (OSA version 6.0). 
Thanks to IBIS large field of view (29° x 29°), SGR 1806-20, 
which is located ~10° from the Galactic Centre, was almost 
constantly observed during the KP. This project was divided in 
two parts: the first one lasted from September 12 to October 
5, 2006, yielding an effective exposure time of ~750 ks on the 
source, while the second one lasted from February 28 to March 
25, 2007 for an exposure of ~550 ks. Our data set consists 
of about 700 individual pointings, but due to the faintness of 
the source we could not extract the spectra from the individual 
pointings. We extracted the images for each pointing in 10 en- 
ergy bands between 20 and 300 keV and added all the individ- 
ual images in order to produce two mosaics (one for each part 
of the KP observation). We then extracted the fluxes from the 
mosaics in each energy band in order to derive the source spec- 
trum, and rebinned the ISGRI response matrix to match our 
energy bands. Since the fluxes were found not to vary within 
the errors between the two parts of the KP, we added the two 
spectra in order to increase our statistics. 

The IBIS/ISGRI spectrum can be well fit Or? = 0.85 
for 5 d.o.f.) with a single power-law, with Y - 1.7 + 0.3, 
and a 20-60 keV flux of (2.8 + 0.4) x lO"" erg cm^^ s '. The 
power-law parameters obtained in the hard X-ray band with 
IBIS/ISGRI and the SuzakulWXT) (Section |Z2T| are consistent 
within the errors. The joint fit of the IBIS/ISGRI and Suzakus 
XIS and HXD-PIN data to the broken power-law plus black- 
body model adopted in Section IZS] yields virtually identical 
best-fit parameters, withx^ = 1-09 for 358 d.o.f.. 

We also fit the IBIS/ISGRI spectrum simultaneously with 
the XMM-Newton/pn spectrum described in the previous sec- 
tion, using a blackbody plus power-law model. The result- 
ing best-fit parameters (F = 1.55 + 0.08, blackbody tempera- 
ture ksT = 0.6 + 0.1 keV and radius Rbb = 2+^ km (at 15 kpc), 
and A^H = 6.6!°'^ x lO^^ cm^^, with;^? = 1-27 for 107 d.o.f.) 



The direction of SGR 1806-20 was scarsely covered by 
INTEGRAL'S X-ray monitor JEM-X jLund et al.l |2003|) . due to its 
smaller (with respect to that of IBIS) field of view (7° diameter). As a 
consequence, only about 200 ks of exposure were available, and this 
was unsufficient for a detection. 




Figure 5. XMM-Newton and INTEGRAL broad-band spectrum 
of SGR 1806-20 (see the online edition of the article for a 
color version of this figure). The data from XMM-Newton/pn 
(black) and INTEGRAL/lBlS (lilac) are fit with the power-law 
(light blu dashed line) plus blackbody (red dot-dashed line) 
model. We also plotted Suzaku/HXD-PIN data using the blu 
circle marks. 



5. Discussion 

The Suzaku, XMM-Newton, and INTEGRAL observations re- 
ported here represent a complementary data set that allows us 
to study the spectral properties of SGR 1806-20 in the broad 
1-100 keV energy range. Although the Suzaku/HXD does not 
have imaging capabilities, we know thanks to INTEGRAL that 
no other bright hard X-ray sources are present in its field of 
view. The uncertainties in the instrumental background (cur- 
rently at the ~5% level) and in the modelling of the Galactic 
Ridge emission are a more relevant concern. Future improve- 
ments in the knowledge of these components may eventu- 
ally allow us to obtain more robust conclusions. Thanks to 
its imaging capabilities, background issues do not affect the 
INTEGRAL observations. However, the IBIS/ISGRI data re- 
quired long integration times, with discontinuous observations 
spanning several months. Thus they provide information on the 
average properties only. The possible presence of long term 
variability was in fact one of the main motivations to per- 
form the Suzaku and XMM-Newton observations simultane- 
ously. With all these caveats in mind we proceed now to dis- 
cuss the broad band spectrum of SGR 1806-20. 

The XMM-Newton and INTEGRAL spectra are consistent 
with an extrapolation of the power-law plus blackbody model 
measured in the 2-10 keV band. Between 12 and 30 keV the 
Suzaku/HXD sensitivity is better than that of INTEGRAL, al- 
lowing to detect SGR 1806-20 during a single 50 ks long ob- 
servation. With respect to die XMM-Newton and INTEGRAL 
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joint fit, the HXD data show an "excess" (see Fig.|5j that can- 
not be completely ascribed to calibration uncertainties between 
the various instruments. 

Given the lack of a direct measure of the Galactic Ridge 
emission around SGR 1806-20, we cannot exclude that this 
excess is due to an underestimation of such contribution to the 
background. If instead the excess is a real feature of the spec- 
trum of SGR 1806-20, its broadband spectrum could be em- 
pirically modeled adopting a power-law with the photon-index 
changing from ~1 to ~2 at ~16 keV, and a blackbody compo- 
nent with kuT ~ 0.8 keV. T his would agree with the results 



reported by' Gotz et al. I (l2006 h. who point out that the hard tails 
of the SGRs are softer than the power-law components mea- 
sured below 10 keV. 

The presence of a down-break in the 10-20 keV spec- 
trum of SGR 1806-20 would have remarkable physical im- 
plications. The soft X-ray emission from magnetar candi- 
dates (SGRs and AXPs) is usually interpreted within the 
twisted magnetosphere model as due to resonant cyclotron up- 
scattering of soft photons from the star surface onto charge s 



flowing into the magnetosphere jThompson et alj 120021) 



Detailed calculations of resonant compto n scattering (RCS) 
spectra have been recently perfor med ( Lyutikov & Gavriij 
2006t Fernandez & Thompson'^2007') and successfully applied 



to fit AXP spectra (Rea et al. 2007). Quit e interestingly, some 
of the model spectra presented by iPernandez & Thompson 
( l2007h exhibit a downward break in the tens of keV range. 
Their overall shape is quite reminiscent of the Suzaku 
XIS/HXD-PIN spectrurn of SG R 1806-20, and, as noted by 
Fernandez & ThompsonI (l2007h . they may also play a role 



in the interpretation of the broadband X-ray spectrum of 
SGR 1900-1-14. In particular, when assuming a (non-thermal) 
top-hat or a broadband velocity distribution for the magneto- 
spheric charges, multiple peaks can appear in the spectrum (see 
their Fig. 6 and Fig. 1 1). The downturn possibly present in our 
data may then be due the presence of a second "hump" (in ad- 
dition to the main thermal one) in the range 10-20 keV. Nobili, 
Turolla, & Zane (in preparation) assuming a 1-D thermal elec- 
tron distribution superimposed to a (constant) bulk velocity, 
found also double humped spectra. In this case the second (and 
only) hump occurs when resonant scattering is efficient enough 
to fill the Wien peak at the temperature of the comptonising 
particles. A spectral break at ~15 keV would translate then in 
a temperature of ~5 keV for the magnetospheric electrons. If a 
more refined treatment of background subtraction confirms the 
spectral break in the X-ray data of SGR 1806-20 this would 
provide important diagnostics for the physical parameters of 
the model. 

In 2003 we started a long-term monitoring program 
to study the time evolution of the spectral properties of 
SGR 1806-20 using the XMM-Newton X-ray satelHte. The 
December 27, 2004 giant flare was a fortunate occuiTence that 
allowed us to observe how the source properties evolved in 
the two years leading up to t he flare and how they changed 
after this drama t ic event (see 



Rea et al.l l2005at iTiengo et al 



Mereghetti et al.1 l2005d. 12007 



2005L for details). The XMM- 



level after the giant flare. A direct comparison of the XMM- 
Newton/pn count rates measured in the different observations 
shows that before the giant flare the flux of SGR 1806-20 in 
the 1-10 keV band was monotonically increasing, while the 
three observations after the flare, and preceding the one re- 
ported here, followed a steady decreasing trend. The September 
2006 observation breaks this long term decay, having a count 
rate higher by 5% with respect to the last XMM-Newton ob- 
servation performed 5 months before. This slight (but statis- 
tically significant) re-brightening might indicate either a tem- 
porary oscillation around an equilibrium flux level or the start 
of a new monotonic flux increase, similar to the one that pre- 
ceded the December 2004 giant flare. This phenomenon was 
interpreted as due to the building up of a magnetospheric twist- 
ing, that determined also the hardening of the X-ray spectrum, 
an inc rease of the spin-down a nd a more intense bursting ac- 
tivity (iMereghetti et al.ll2005cl) . The relatively high burst rate 
observed during the XMM-Newton and Suzaku observations of 
September 2006 (see Section lZST i is therefore another indica- 
tion of a possible increase of the magnetospheric twisting in 
SGR 1806-20, but, before a new XMM-Newton observation 
will be performed, only the monitoring of the frequency and 
intensity of SGR 1806-20 bursts can tell us if the evolution is 
erratic or follows a stable trend. The recent report of a bright 
burst from SGR 1806-20 iGolenetskii et al.ll2007l;|Perotti et al " 



2007h seems actually to favour the second hypothesis. 

Two of the bursts detected during the Suzaku observation 
were bright enough to allow spectral analysis. In both cases, the 
broadband spectrum (2-100 keV) revealed the presence of two 
components: a soft component which is well reproduced by a 
blackbody with ksT ~ 2-4 keV, and a harder one whose spec- 
tral shape is not firmly established and can be equally well fit 
with a power-law, a hot bremsstrahlung or a second blackbody 
(See Table|2]and Fig.|4]l. In absence of robust theoretical pre- 
dictions, we can not exclude that a two component model sim- 
ply reflects our ignorance of the correct spectral shape, and has 
therefore a purely phenomenological significance. However, it 
is worth noticing that, from their recent analysis of a sample of 
50 bursts detected from SGR 1806-20 with HETE-2 from 2001 



to 2005. lNakagawa et alJ (l2007h have suggested the presence of 
a time delay between the 30-100 keV and the 2-10 keV emis- 
sion. Although such a delay can be attributed to an intrinsic, 
rapid spectral softening, an alternative, and simpler, interpreta- 
tion invokes the presence of two separate emitting regions. 

Let us consider a scenario in which the two components 
are physically distinct and let us consider the hard compo- 
nent first. In the magnetar scenario, short bursts are usually as- 
cribed to either reconnection phenomena in the external mag- 
netosphere (eve ntually modulated by a tearing instability, see 
Lyutikovll2003h or movements of the footprints of the exter- 
nal magnetic field, produced by crustal deformations or frac- 
ture s driven by the stress exerted by the internal field helic - 
ity ([Thompson & Duncan|[l99l bOOlL iTTTompson et al.ll2002l) . 



Newton data showed a doubling of the flux in the September- 
October 2004 followed by a gradual recovery to the "historical" 



Both kind of processes lead to the generation and launch of 
an Alfven wave, which produces and accelerates a particle 
cascade, and ultimately is detected as a burst. The emerging 
spectrum is expected to be synchrotron dominated, unless the 
Alfven wave is temporarily trapped in a fireball and thermal- 
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ized. Therefore, both the BB + PL and BB + BB spectral fits 
are consistent with a scenario in which such an Alfven wave is 
responsible for the hard component. We notice that, although 
a fireball formation is not required to explain short bursts (and 
therefore usually not invoked in such cases), to our knowledge 
there is no a priori reasons why a small fireball can not be cre- 
ated and evaporated in a sub-second time scale, giving rise to 
a thermal spectrum. A point in favour of this interpretation is 
that, in the BB + BB fit, the temperature of the hot black- 
body is remarkably close to the minimal temperature above 
which a fireball thermalizes, kuT ^ 1 1 (j ^Ns/lOkm)''^^ keV, 



according t o Thomps on & DuncanI (12001 ). eq. [71] (see also 



OUve et al 



2004, for similar findings based on longer dura- 
tion burst). The third model of the hard component which 
is compatible with our data is a bremsstrahlu ng emission a t 
~10 keV. Quite recently, [T hompson & Beloborodov (l2005b 
and iBeloborodov & ThompsonI (l2007l) discussed the electro- 
dynamics of the magnetar coronae and the production mech- 
anisms for soft gamma-rays. In particular, their model pre- 
dicts the existence of a thin transition layer between the corona 
and the thermal photosphere, where Langmuir turbulence can 
be excited by a downward beam of current-carrying charges. 
As a result, the transition layer can be heated up to a typical 
temperature of -100 keV, and emit, approximatively, an op- 
tically thin bremsstrahlung at a s ingle temperature. Although 
Thompson & Beloborodovl (l2005a model was originally devel- 
oped in connection with the persistent hard emission of magne- 
tars, the predicted bremsstrahlung temperatures are remarkably 
close to those we detected during the two bursts, suggesting 
that a similar mechanism may instead be activated during peri- 
ods of activity. 

Our results about the spectral modelling of the soft X-ray 
component are more robust, inasmuch as all our spectral fits 
require the presence of a cold blackbody with ~ 2-4 keV. 



This i s in agreement with s imilar findings by lOlive et al . 
( |2004 . iNakagawa et al.1 (|2007|) . and lFeroci et"aLl (|2004 . This 



component is usually interpreted as due to emission from a 
fraction of the star surface (which can be as large as the whole 
star in the case of our BB + BB fit) heated by returning cur- 
rents. Alternatively, it has been suggested that the soft compo- 
nent may originate up in the magnetosphere (<700 km), pre- 
suma bly due to a delayed emission process d Nakagawa et al. 
20071). Here we only notice that, although the spectra of our 



two events are compatible with emission from the star sur- 
face, the radius of the cold blackbody as measured during other 
short bursts can reach values much higher than 50-100 km 
( Nakagawa et al J|2007 , similar findings have been found in the 
case of SGR 1900-H14 bursts measured with Swift, Israel et al. 
(in preparation)). One possible explanation is that part of the 
flare energy is intercepted and reprocessed in a larger region 
and re-emitted at a lower temperature. In such scenario, the ra- 
dius of the reprocessing region can then vary depending on the 
fraction of material that is i ntercepted, and is not bounded by 
the value of the star radius. Thompson & DuncanI (12001,) con- 
sidered the equilibrium state of a pair corona sufficiently ex- 
tended that the local value of the magnetic field is B «: Bqed, 
so that photon splitting can be ignored (if the magnetic field 
scales as a dipole, this occurs above ~3 star radii, for a polar 



surface value of ~10'^ G). They found that a stable balance 
between heating and diffusive cooling requires a continuous 
source of ordinary photons that can be provided, for instance, 
by external illumination. If the corona intercepts part of the 
flare beam (which in their treatment is assumed to originate in 
a trapped fireball, although in general this is not necessarily 
required), equihbrium is possible below a critical luminosity 
given by 



L<L„. = 1.5xlO«.-|2^ 



1 10 km/ 



erg s 



where Tg is the pair temperature in the corona, teo ~ 1 is the 
scattering depth for ordinary to extraordinary mode conversion 
and R the radius of th e emitting part of the corona (see equa- 
tions [84] and [89] in Thompson & Duncan (1200 1") and note 
that there is a typo in their equation [89]: it should contain 
instead of R^^). Even for an emission region as small as 
5-20 km (as inferred by our best fit of the low temperature 
blackbody) and a temperature of kgT ~ 2-4 keV, this is well 
above the luminosity emitted during the two bursts detected by 
Suzaku. Therefore, simply on the basis of energetics, relatively 
large emitting regions for the cold blackbody are compatible 
with the temporary formation of a pair corona, sustained by a 
fraction of the flare energy. When the heating rate ceases, the 
pair atmosphere contracts and quickly evaporates. In order to 
derive firmer conclusions a more detailed analysis is needed, 
mainly in assessing the possibility that the intercepted beam 
is thermalized and re-emitted as a blackbody. This study is be- 
yond the purpose of this paper, and will be presented elsewhere 
(Israel et al. 2007, in preparation). 

Acknowledgements. This work is based on data from observations 
with Suzaku, INTEGRAL, and XMM-Newton. Suzaku is a Japan's 
mission developed at the Institute of Space and Astronautical Science 
of Japan Aerospace Exploration Agency in collaboration with U.S. 
(NASA/GSFC, MIT) and Japanese institutions. We would like to 
thank the Suzaku operations team for their support. XMM-Newton 
is an ESA science mission with instruments and contributions di- 
rectly funded by ESA Member States and NASA. INTEGRAL is 
an ESA project with instruments and science data centre funded by 
ESA member states (especially the PI countries: Denmark, France, 
Germany, Italy, Switzerland, Spain), Czech Republic and Poland, and 
with the participation of Russia and the USA. This research has made 
use of HEASARC online services, supported by NASA/GSFC, and 
of the SIMBAD database, operated at CDS, Strasbourg, France. The 
Italian authors acknowledge the support of the Italian Space Agency 
(contracts 1/023/05/0 and 1/008/07/0) and the Italian Ministry for 
University and Research (grant PRIN 2005 02 5417). SZ acknowl- 
edges STFC for support through an Advanced Fellowship. DG ac- 
knowledges the French Space Agency (CNES) for financial support. 
NR is supported by an NWO Post-doctoral Fellowship. Finally, we 
wish to thank the anonymous referee for providing useful comments 
on the manuscript. 

References 

Anders, E. & Grevesse, N. 1989, Geochim. Cosmochim. Acta, 
53, 197 

Bamba, A., Koyama, K., Hiraga, J. S., et al. 2007, PASJ, 59, 
209 



p. Esposito et al.: X-ray observations of SGR 1806-20 



11 



Beloborodov, A. M. & Thompson, C. 2007, ApJ, 657, 967 
Cameron, P. B., Chandra, P., Ray, A., et al. 2005, Nature, 434, 
1112 

Corbel, S., Wallyn, P, Dame, T. M., et al. 1997, ApJ, 478, 624 
Fenimore, E. E., Laros, J. G., & Ulmer, A. 1994, ApJ, 432, 742 
Fernandez, R. & Thompson, C. 2007, ApJ, 660, 615 
Feroci, M., Caliandro, G. A., Massaro, E., Mereghetti, S., & 

Woods, R M. 2004, ApJ, 612, 408 
Frederiks, D. D., Golenetskii, S. V., Palshin, V. D., et al. 2007, 

Astronomy Letters, 33,1 
Gaensler, B. M., KouveUotou, C, Gelfand, J. D., et al. 2005, 

Nature, 434, 1104 
Golenetskii, S., Aptekar, R., Mazets, E., et al. 2007, GCN Circ. 

6684 

Gotz, D., Mereghetti, S., Tiengo, A., & Esposito, P. 2006, 

A&A,449, L31 
Gruber, D. E., Matteson, J. L., Peterson, L. E., & Jung, G. V. 

1999, ApJ, 520, 124 
Hurley, K., Boggs, S. E., Smith, D. M., et al. 2005, Nature, 434, 

1098 

Hurley, K., Cline, T, Mazets, E., et al. 1999, Nature, 397, 41 
Israel, G., Covino, S., Mignani, R., et al. 2005a, A&A, 438, LI 
Israel, G. L., Belloni, T, Stella, L., et al. 2005b, ApJ, 628, L53 
Kokubun, M., Makishima, K., Takahashi, T, et al. 2007, PASJ, 
59,53 

Kosugi, G., Ogasawara, R., & Terada, H. 2005, ApJ, 623, L125 
Kouvehotou, C, Dieters, S., Strohmayer, T., et al. 1998, 

Nature, 393, 235 
Koyama, K., Tsunemi, H., Dotani, T, et al. 2007, PASJ, 59, 23 
Laros, J. G., Fenimore, E. E., Fikani, M. M., Klebesadel, R. W., 

& Barat, C. 1986, Nature, 322, 152 
Laros, J. G., Fenimore, E. E., Klebesadel, R. W., et al. 1987, 

ApJ, 320, LI 11 

Lebrun, R, Leray, J. P, Lavocat, P., et al. 2003, A&A, 411, 
L141 

Lebrun, R, Terrier, R., Bazzano, A., et al. 2004, Nature, 428, 
293 

Lund, N., Budtz-j0rgensen, C, Westergaard, N. J., et al. 2003, 

A&A,411,L231 
Lyutikov, M. 2003, MNRAS, 346, 540 
Lyutikov, M. & Gavriil, R R 2006, MNRAS, 368, 690 
Mazets, E. P., Golentskii, S. V., Ilinskii, V. N., Aptekar, R. L., 

& Guryan, I. A. 1979, Nature, 282, 587 
McClure-Griffiths, N. M. & Gaensler, B. M. 2005, ApJ, 630, 

L161 

Mereghetti, S., Esposito, P, & Tiengo, A. 2007, Ap&SS, 308, 

13 

Mereghetti, S., Gotz, D., Mirabel, I. R, & Hurley, K. 2005a, 
A&A, 433, L9 

Mereghetti, S., Gotz, D., von Kienlin, A., et al. 2005b, ApJ, 
624, L105 

Mereghetti, S., Tiengo, A., Esposito, P., et al. 2005c, ApJ, 628, 

938 

Mitsuda, K., Bautz, M., Inoue, H., et al. 2007, PASJ, 59, 1 
Murakami, T., Tanaka, Y., Kulkami, S. R., et al. 1994, Nature, 

368, 127 

Nakagawa, Y. E., Yoshida, A., Hurley, K., et al. 2007, PASJ, 
59, 653 



Olive, J.-R, Hurley, K., Sakamoto, T, et al. 2004, ApJ, 616, 
1148 

Perotti, R, Giuhani, A., Mereghetti, S., et al. 2007, GCN Circ. 
6688 

Rea, N., Israel, G., Covino, S., et al. 2005a, Astron. Tel. 645 
Rea, N., Tiengo, A., Mereghetti, S., et al. 2005b, ApJ, 627, 
L133 

Rea, N., Zane, S., Lyutikov, M., & Turolla, R. 2007, Ap&SS, 
308, 61 

Serlemitsos, P J., Soong, Y, Chan, K.-W., et al. 2007, PASJ, 

59,9 

Striider, L., Briel, U., Dennerl, K., et al. 2001, A&A, 365, L18 
Sugizaki, M., Mitsuda, K., Kaneda, H., et al. 2001, ApJS, 134, 
77 

Takahashi, T, Abe, K., Endo, M., et al. 2007, PASJ, 59, 35 
Taylor, G. B., Gelfand, J. D., Gaensler, B. M., et al. 2005, ApJ, 

634, L93 

Terasawa, T, Tanaka, Y. T, Takei, Y, et al. 2005, Nature, 434, 
1110 

Thompson, C. & Beloborodov, A. M. 2005, ApJ, 634, 565 
Thompson, C. & Duncan, R. C. 1995, MNRAS, 275, 255 
Thompson, C. & Duncan, R. C. 2001, ApJ, 561, 980 
Thompson, C, Lyutikov, M., & Kulkami, S. R. 2002, ApJ, 574, 

332 

Tiengo, A., Esposito, P., Mereghetti, S., et al. 2005, A&A, 440, 
L63 

Turner, M. J. L., Abbey, A., Arnaud, M., et al. 2001, A&A, 
365, L27 

Ubertini, P, Lebrun, P., Di Cocco, G., et al. 2003, A&A, 411, 
L131 

Valinia, A. & Marshall, R E. 1998, ApJ, 505, 134 
Woods, R M., KouveUotou, C, Ringer, M. H., et al. 2007, ApJ, 
654, 470 



